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Abstract 



• We compute nuclear spin dependent structure functions using a dynamical model for 



bound nucleon densities and hence calculate nuclear modifications to asymmetries ob- 
served in recent doubly polarised deep inelastic scattering experiments. We conclude 



' that while the individual densities are changed substantially by nuclear effects, the 

asymmetries themselves are largely insensitive to these changes. 
T— I ■ 

■ Recently a model was proposed [|I| to explain the observed differences between free nu- 

■ cleon and bound nucleon structure functions in deep inelastic lepton nucleon scattering 
^ . (DIS). This model used a dynamical approach, involving modifying a free nucleon input 
O '• density distribution at a low input scale, = fi"^ = 0.23 GeV^, due to nuclear effects, and 
^ I then evolving the resultant modified bound nucleon densities to the required Q"^ scale of the 

■ experiment. The model gave satisfactory agreement with available data in a fairly broad 
p^- range, from 0.5-30 GeV^. 

p . It is interesting to ask how this model can be extended to a study of spin dependent 

D ! bound-nucleon densities. The question is not merely academic as, in fact, data on the spin 
^ ' dependent deuteron and neutron structure functions have been obtained [0 from deuteron 
and ^He targets. Nuclear effects in deuteron are known to be small (though measureable) , 
^ ■ since the deuteron is a loosely bound nucleus. There have been a number of papers |§ 
dealing with nuclear modifications of spin asymmetries and structure functions in the case 
of the deuteron^. We therefore confine our attention to possible nuclear effects on the double 
spin asymmetry measurements made with helium nuclei. In this case, it was pointed out 
by Woloshyn 0| that the protonic contribution to the asymmetry is negligible so that the 
^He double spin asymmetry is sensitive to the spin dependent neutron structure function, 
gi{x, Q^). However, there may be additional modifications due to the presence of the nuclear 
medium, which we propose to study here. These are especially of importance for checking 
the validity of the Bjorken Sum rule. Our main conclusion is that the individual (spin 
independent as well as spin dependent) structure functions undergo substantial modifications 
due to nuclear effects; however, their ratio — the asymmetry — which is the measured quantity, 
is largely free from these and so gives hope that the neutron structure function may be 
unambiguously determined from such a measurement. 

^ e-mail : indu@hal 1 . physik. uni-dortmund . de 

^Depending on the model, corrections due to nuclear effects in deuterium can be as large as 10%. 
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The quantity of interest is the double spin asymmetry, 



where gi and Fi are the spin dependent and spin independent structure functions corre- 
sponding to the nucleus A. We are therefore interested in studying possible deviations of 
the measured asymmetry, Af*^, from the required neutron asymmetry, A", due to nuclear 
effects. 

As a starting point we note that the corresponding study of unpolarised bound nucleon 
densities used the GRV density parametrisation as an input. We shall therefore use the 
GRVs spin dependent densities as an input in the corresponding polarised problem. This 
is essential if we are to retain the definition, 

qf{x) = qj{x) + qj{x) ; ^2) 
q^(x) = qf{x) - qj{x) , 

where q^ and qj are the positive and negative helicity densities of /-flavour quarks (and 
similarly for gluons) in either free or bound nucleons. 



1 The Spin Independent Nuclear Densities 

We now quickly review the model before we apply it to the polarised case. This is also useful 
as the polarised case essentially follows along the lines of the unpolarised problem. The free 
nucleon input densities are modified by both nuclear swelling and binding effects. Nucleon 
swelling causes not only a depletion of all parton densities at large- and small-x, but also an 
enhancement at intermediate-x 0. The relative increase in the nucleon's radius is 6a, where 
{Rn + AR{A))/Rn = 1 + 6a, and is given by, 

6a = [I- Ps{A)]6,oi + Ps{A)6,oi/2 . 

The second term corrects for surface effects in the usual manner. Here Ps{A) is the proba- 
bility of finding a nucleon on the nuclear surface while 6^1 parametrises the swelling of the 
nucleon in the interior of a heavy nucleus and is the only free parameter in the calculation. 
It was fixed to be 5^1 = 0.15 in the unpolarised calculation [Q. 

The distortions of the density distributions due to swelling, being purely geometrical, 
conserve the total parton number and momentum of each parton species (i.e., the first 
and second moments of the distributions are unchanged). Furthermore, the third moments 
are modified in a well- determined way. Specifically, the first three moments of the parton 
distributions in a free (^at) and bound (g^) nucleon at /i^ are related by 

(?a(/u2))i = (giv(/u2))i , 

(9a(^^))2 = (g7v(/i^))2 , /ox 

{{<!N{f^'))s-{qNif^'))lY^' _ 



((gA(/))3-(gA(^^)) 



2N\2Nl/2 
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The first two equations imply number and momentum conservation of partons and tlie last 
incorporates the swelling effect 0. 

These are then used as constraint equations to determine the bound-nucleon densities in 
terms of the known free-nucleon ones. We use the Gliick, Reya and Vogt (GRV) parametri- 
sation |Q of the free-nucleon distributions at = fi^ = 0.23 GeV^. We find these most 
appropriate for our purpose as each of their input densities is integrable (there are finite 
number of partons at /i^). 

It is now possible to determine the bound nucleon densities (for q = valence quarks, 
uv, dv, sea quarks, S, and gluons, g) in terms of the free densities, qN- We parametrise the 
free as well as bound nucleon distributions in the form, 

q{x) = Nx'^il-xfPix), (4) 

where P{x) is a polynomial. We take PA,q{x) = PN,q{x), for simplicity. Then the changes 
in the three main parameters, N, a, and (3, due to swelling, and hence the bound-nucleon 
densities, are immediately determined by the constraints in eq (3). This fixes the input 
bound-nucleon densities at = fi^. 

We now discuss the binding effect. The attractive potential describing the nuclear force 
arises from the exchange of mesons. Hence the energy required for binding is taken away 
solely from the mesonic component of the nucleon, and not from its other components. At 
the starting scale, = /x^, we identify these mesons to be just the sea quarks in the bound 
nucleon. Hence, the momentum fraction carried by the sea quarks in a nucleon bound in a 
nucleus at = /i^ will be reduced. The extent of reduction is determined by the binding 
energy per nucleon [|l|, which is given by the well-known Weizacker mass formula. 

These nuclear effects completely determine the spin independent input bound-nucleon 
densities. These are then evolved using the usual Altarelli Parisi evolution equations to 
obtain the densities at a required value of Q^. 

At the time of interaction, there is a further depletion of the sea densities, which occurs 
whenever there is nucleon-nucleon interaction, caused by parton-nucleon overlap. When a 
parton having a momentum fraction, x, of the parent nucleon momentum, P^, is struck, it 
is off-shell and localised to a distance, AZ ~ l/(2xP/v) (in the Breit frame). For sufficiently 
small X, AZ becomes large and can exceed the average 2-nucleon separation^. 

The struck parton must return to the parent nucleon within the interaction time, as 
required by the uncertainty principle. However, while it extends outside the parent nucleon, 
it can interact with other nucleons in the nucleus. Such an interaction between two nucleons 
caused by parton-nucleon overlap results in loss of energy of the parent nucleon, mimicking 
exactly the effect of binding. Hence we call this the second binding effect and assume its 
strength to be the same as that due to the usual binding. This immediately fixes the loss 
in sea quarks (valence quarks are not depleted due to the requirement of quantum number 
conservation) due to this effect to be 

S'^{x,Q^) = K'{A)Sa{x,Q^) , (5a) 

•^Although the spatial extent of a single coloured parton cannot exceed the range of QCD confinement, the 
struck parton can combine with a wee parton and form a colourless scalar with vacuum quantum numbers 
which can then escape from the nucleon. 
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where the depletion factor is, 



K\A) 



1, 

1 — 2I3{xqx^ 
1 - 2(3{xqXa 



when X > Xq] 
1), when xa < X < Xq] 
1), when X < xa, 



1 



(56) 
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where /3 is the same as in usual binding, viz.. 



0.037/2 , 



(6) 



Mjv(^jv(/i2))2 



U{^^) being the binding energy between each pair of nucleons, which is known. The lim- 
iting values, Xo = l/(2M7v(iAr) (where dj^ is the average correlation distance between two 
neighbouring nucleons in the lab frame), and xa = ^/ {^RaMn) (where 2Ra — IARa is the 
average thickness of the nucleus), determine the starting and saturation values respectively 
of this shadowing effect; the latter occurs when the struck quark wave function completely 
overlaps the nucleus in the 2;-direction. In general, a parton with a momentum fraction, x, 
Xa ^ X < Xo, can overlap {n — 1) other nucleons, where n = 1/ {2Mj^d]yx) = Xq/x. Due 
to the applicability of the superposition principle to the scalar field interaction with various 
nucleons, the loss of energy due to interaction with each of the nucleons over which the struck 
quark wave function extends, is equal and additive, and thus explains the depletion factor 
in eq (5). Since this effect acts on the intermediate state of the probe-target interaction, it 
does not participate in the QCD evolution of the initial state. 

Nuclear modification due to binding and swelling at the input scale = fi"^, and parton- 
nucleon overlap due to the second binding effect at the scale of the scattering together 
determine the structure function, Ff{x, Q"^), of a nucleon bound in a nucleus A. The model 
gives good agreement with available data 

We now proceeed to an analysis of the corresponding spin dependent densities. 

2 The Spin Dependent Nuclear Densities 

The same nuclear effects of binding and swelling affect the spin dependent densities also. 
This is because they influence the positive and negative helicity densities, out of which the 
spin independent and spin dependent densities are composed (see eq (2)). The entire swelling 
effect can now be rephrased as the effect of swelling on individual helicity densities, so that 
equations analogous to (3) are valid for the spin dependent densities, q{x), as well. This can 
be seen as follows: Swelling simply rearranges the parton distributions in the bound nucleon; 
there is no change in the number of each parton species. In particular, each helicity type is 
also conserved, i.e.. 



Hence, their sum and difference is also conserved. The former is contained in the first 
equation of the equation set (3); the latter implies, for the polarised combination. 
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Note that {q{li^))n = (?^(yU^))n + (?~(/^^))n for every moment, n, for both the free and 
bound nucleon, and similarly for the spin dependent density as well. Similarly, since the 
momentum carried by each helicity density is unchanged, momentum conservation between 
the free and bound nucleon also holds for the sum and difference of the helicity densities. 
The corresponding equation for the sum is the second equation in (3); the equation for the 
helicity difference is 

{ui^^'))2 = mii'))2 . {7b) 

The extension of the third of the equations in (3) to the spin dependent case is not as 
straightforward. Every helicity density, qf{x), {h = +,— ), spreads out over a larger size, 
or, equivalently, gets pinched in momentum space, according to Heisenberg's uncertainty 
relation, ApAx = 1. Applying this to each helicity type, for each flavour, we have, 

((4(/i^))3-(4(/i^))2)^^^ , , iJMll)h^JMllM}ll = ^ ^ X 

However, for later convenience, we prefer to use analogous expressions for the sum and 
difference, qf and g/, rather than for the individual helicity densities. Hence, the third of the 
constraints arising from swelling, i.e., the third of eq (3) and its spin dependent counterpart 
read, 



1 + Sa. (7c) 



The error involved between the exact expressions, eq (8), and their approximations, eq (7c), 
is a term proportional to (1 — (1 + Sa)'^) and is of order 6a- This term mixes spin dependent 
and spin independent moments; however, since 6a is small (about 10%), these errors are 
small, and can be ignored. We are therefore justified in using eq (7c) rather than eq (8) to 
constrain the second moments of the parton densities. The three sets of equations, (7a-c), 
thus provide the three sets of constraint equations, analogous to the set (3), with which we 
can fix the input bound nucleon spin dependent densities. 

The modified input densities are thus determined, given a set of valid input free nucleon 
distributions, which we take to be the Gliick, Reya, and Vogelsang 'standard' set (GRVs) 
0. These densities can also be parametrised in a form similar to eq (4); in fact, every 
spin dependent density is a factor of the form of the RHS of eq (4) times the corresponding 
unpolarised density. Hence there are again three constraint equations which serve to fix the 
three main parameters, a, (3, and for the corresponding bound nucleon spin dependent 
densities. 

Binding causes loss of energy in the sea: this is due to loss of mesons from the nucleon. 
Since these mesons are spin-0 bosons, it is clear that no spin is lost from the sea due to 
binding (equal numbers of positive and negative helicity partners are lost). Hence we see 
that binding changes the sum, but not the difference of the helicity densities^. 

We thus obtain the input polarised densities analogous to the unpolarised ones. These 
are then evolved to the scale of interest. 

"^It is possible that p, etc., mesons also participate in this interaction, leading to a change in the polarised 
sea densities, but this component is small and we neglect it. 



5 



At the time of interaction, the second binding effect apphes to struck partons with mo- 
mentum fraction x < Xq, as in the unpolarised case. The mechanism for this depletion 
is independent of the hehcity of the quark, and so this effect is identical in both the spin 
independent as well as spin dependent cases. Hence, the spin dependent structure function, 
gi{x, Q^) can now be computed by evolving the modified input spin dependent densities to 
the required value of Q^, and including the second binding effect. 

Finally, we display the equivalent neutron bound-nucleon structure functions at an arbi- 
trary scale, > jj? (with R = ai/ax = 0) : 



"M/'^ n2\ _ 1 



Fr{x,Q 



g'l'\x,Q') = ^ 



u^{x, g2) + 4 d^{x, + k'{a)Sa{x, q 
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(9) 



q^{x,fi'^) incorporates the effect of swelling on every input parton density, q^^Xyfi"^), as 
well as that of binding for the unpolarised densities, and the corresponding Q^-dependent 
quantities that appear here are these input densities, evolved suitably to the required scale. 
K'{A) incorporates the second binding effect, at Q^, as discussed above. The experimentally 
measured asymmetry, and quantity of interest, are the ratios, at the scale Q^, for the neutron 
bound in the helium nucleus and for a free neutron: 

n/He n 
/I meas 9l /I rend 9l 

- n/He ' ~ pn ' l-^UJ 

and can thus be computed. (We use the free and bound nucleon unpolarised structure 
functions from [Q). Note that the input spin dependent densities (which are taken from 
0) were actually fitted to both the free proton as well as deuteron and ^He spin dependent 
data; however, we use them here as the free nucleon parametrisations (which is permissible 
especially in view of the large error bars on presently available data). Furthermore, the 
smearing effect of Fermi motion (at large x) is neglected in this work for simplicity. Hence 
our results are not valid at large x. 

In fig. 1 we give the results of our computations for the measured (bound nucleon) 
and required (free nucleon) spin dependent structure function, for typical values of Q^, 
Q"^ = 1,4 GeV^. We see that the deviations of the bound neutron structure function can be 
as large as 10-15% at small x and about 6% at intermediate x values. The data points plotted 
on this graph correspond to the values extracted at = 4 GeV^ from a measurement of the 
asymmetry by the E142 Collaboration [Q] (with R = 0) and indicate the size of the error bars 
in currently available data. In fig. 2, we plot the asymmetries at = 4 GeV^. The data 
points here correspond exactly to the E142 data and therefore go over a range of with a 
mean of about 2 GeV^; however, the asymmetry is not very sensitive to in the x range of 
the available data. Notice that in this case, nuclear effects cause not more than 5% deviation 
in the asymmetry at both small and intermediate values of x. The deviation is slightly larger 
at larger x, x > 0.4, but this is due to the fact that the neutron spin dependent structure 
function changes sign near this value, and hence this deviation cannot be considered to be 
significant. 

In short, we see that nuclear effects, though significant, equally affect both the spin de- 
pendent as well as the spin independent structure functions in such a way that the measured 
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asymmetries are to a great extent independent of them. Since it is the asymmetry rather 
than the structure function which is measured in a polarised experiment, much smaller errors 
on data are required before these small deviations due to nuclear effects become observable 
in such experiments. On the other hand, as already stated, this seems to make possible 
clean and unambiguous extraction of the relevant free nucleon structure functions from a 
measurement of double spin asymmetries with such light nuclear targets. 

Acknowledgements: I thank M. Gliick for suggesting the idea that a computation of 
nuclear effects on spin structure function data may be interesting; I thank E. Reya for a 
critical reading of the manuscript, and M. Stratmann for providing the relevant FORTRAN 
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Figure Captions 

Fig. 1 The free and bound nucleon spin dependent structure function for = 1, 4 GeV^ as 
a function of x are shown as solid and dashed lines respectively. The structure function data 
are extracted at = 4 GeV^ from the asymmetries measured by the E142 collaboration. 

Fig. 2 The bound and free nucleon asymmetries for Q'^ — A GeV^ as a function of x are 
shown as solid and dashed lines respectively. The data are from the E142 collaboration. 
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Fig. 1 
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